Geographical networks stochastically constructed 
by a self-similar tiling according to population 
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In real communication and transportation networks, the geographical positions of nodes are very 
important for the efficiency and the tolerance of connectivity. Considering spatially inhomogeneous 
positions of nodes according to a population, we introduce a multi-scale quartered (MSQ) network 
that is stochastically constructed by recursive subdivision of polygonal faces as a self-similar tiling. 
It has several advantages: the robustness of connectivity, the bounded short path lengths, and the 
shortest distance routing algorithm in a distributive manner. Furthermore, we show that the MSQ 
network is more efficient with shorter link lengths and more suitable with lower load for avoiding 
traffic congestion than other geographical networks which have various topologies ranging from river 
to scale-free networks. These results will be useful for providing an insight into the future design of 
ad hoc network infrastructures. 

PACS numbers: 89.75.Hc, 02.50.Ga, 89.20.Ff, 89.40.-a 



I. INTRODUCTION 

Since the common topological characteristics called 
small-world (SW) and scale-free (SF) have been re- 
vealed in social, technological, and biological networks, 
researches on complex networks have attracted much at- 
tention in this decade through the historical progress [ ] . 
The topological structure is quite different from the con- 
ventionally assumed regular and random graphs, and has 
both desirable and undesirable properties: short paths 
counted by hops between any two nodes and the fault- 
tolerance for random node removals [2, 3] on the one 
hand but the vulnerability against intentional attacks on 
hubs [3, 4] on the other hand. In the definition of a net- 
work, spatial positions of nodes and distances of links 
are usually ignored, these simplifications arc reasonable 
in some networks such as the World- Wide- Web, citation 
networks, and biological metabolic networks. However, 
as real-life infrastructures, in communication networks, 
transportation systems, and the power grid, they are cru- 
cial factors; a node is embedded in a mixing of sparse and 
dense areas according to the population densities, and 
a connection between nodes depends on communication 
efficiency or economical cost. Thus, a modeling of geo- 
graphical networks is important to understand the fun- 
damental mechanism for generating both topological and 
spatial properties in realistic communication and trans- 
portation systems. 

Many methods for geographically constructing com- 
plex networks have been proposed from the viewpoints 
of the generation mechanism and the optimization. As 
a typical generation mechanism, a spatially preferential 
attachment is applied in some extensions [5-9] of the 
Barabasi- Albert (BA) model [10]. As typical optimiza- 
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tions in the deterministic models [11—14], there are sev- 
eral criteria for maximizing the traffic under a constraint 
[11], minimizing a fraction of the distance and node de- 
gree with an expectation of short hops [12], and mini- 
mizing a sum of weighted link lengths w.r.t the edge be- 
tweenness as the throughput [] 3] or the forwarding load 
at nodes [14]. In these methods, various topologies rang- 
ing from a river network to a SF network emerge accord- 
ing to the parameter values. The river network resembles 
a proximity graph known in computer science, which has 
connections in a particular neighbor relationship between 
nodes embedded in a plane. In constructing these net- 
works, it is usually assumed that the positions of nodes 
are distributed uniformly at random, and that a pop- 
ulation density or the number of passengers is ignored 
in communication or transportation networks except in 
some works [13-15]. However, in real data [16], a popula- 
tion density is mapped to the number of router nodes on 
Earth, the spatial distribution of nodes does not follow 
uniformly distributed random positions represented by a 
Poisson point process. Such a spatially inhomogeneous 
distribution of nodes is found in air transportation net- 
works [17] and in mobile communication networks [18]. 

On the other hand, geometric methods have also been 
proposed as another generation mechanism, in which 
both SW and SW structures are generated by a recursive 
growing rule for the division of a chosen triangle [19-22] 
or for the attachment aiming at a chosen edge [23-25] in 
random or hierarchical selection. The position of a newly 
added node is basically free as far as the geometric oper- 
ations are possible, and has no relation to a population. 
Thus, the spatial structure with geographical constraints 
on nodes and links has not been investigating enough. In 
particular, considering the effects of a population on a ge- 
ographical network is necessary to self-organize a spatial 
distribution of nodes that is suitable for socio-economic 
communication and transportation requests. 

In this paper, as a possibility, we pay attention to a 



combination of complex network science and computer 
science (in particular, computational geometry and rout- 
ing algorithm) approaches. This provides a new direc- 
tion of research on self-organized networks by taking into 
account geographical densities of nodes and population. 
We consider an evolutionary network with a spatially in- 
homogeneous distribution of nodes based on a stochastic 
point process. Our point process differs from the tessella- 
tions for a Voronoi partitioning with different intensities 
of points [2G] and for a modeling of crack patterns [27]. 
We aim to develop a future design method of ad hoc net- 
works, e.g., on a dynamic environment which consists of 
mobile users, for increasing communication requests, and 
wide-area wireless and wired connections. More precisely, 
the territory of a node defined as the nearest access point 
is iterativcly divided for load balancing of communication 
requests which are proportional to a population density 
in the area. A geographical network consisting of a self- 
similar tiling is constructed by recursive subdivision of 
faces according to a population. It is worth noting that 
positions of nodes and a network topology are simulta- 
neously decided by the point process in a self-organized 
manner. Furthermore, the geographical network has sev- 
eral advantages [28] : the robustness of connectivity, the 
short path lengths, and the decentralized routing algo- 
rithm [29]. Taking these advantages into consideration, 
we generalize the point process biased by a population 
for constructing a geographical network, and investigate 
the traffic load on the shortest distance routing. 

The organization of this paper is as follows. In Sec. 
II, we introduce a more general network model for self- 
similar tilings than the previous model [28] based on tri- 
angulations. By applying the geometric divisions, we 
construct a geographical network according to a given 
population. In Sec. Ill, we show the properties of the 
shortest path and the decentralized routing without a 
global table for packet transfers as applied in the Inter- 
net. We numerically investigate the traffic load in the 
proposed network, comparing it with the load in other 
geographical networks. In particular, we show that our 
geographical network is better than the state-of-the-art 
geographical networks in terms of shorter paths and link 
lengths, and of lower load for avoiding traffic congestion. 
In Sec. IV, we summarize the results and briefly discuss 
further studies. 



II. MSQ NETWORK MODEL 

We introduce a multi-scale quartered (MSQ) network, 
which is stochastically constructed by a self-similar tiling 
according to a given population. Let us consider the basic 
process of network construction [28]. Each node corre- 
sponds to a base station for transferring packets, and a 
link between two nodes corresponds to a wireless or wired 
communication line. Until a network size N is reached, 
the following process is repeated from an initial configu- 
ration which consists of equilateral triangles or squares. 



At each time step, a triangle (or square) face is chosen 
with a probability proportional to the population in the 
space of the triangle (or square). Then, as shown in 
Figs. 1(a) and 1(b), four smaller triangle (or square) 
faces are created by adding facility nodes at the inter- 
mediate points on the communication links of the chosen 
triangle (or square). This process can be implemented 
autonomously for a division of the area with the increase 
of communication requests. Thus, a planar network is 
self-organized on a geographical space. Figure 2 shows 
an example of the geographical MSQ network according 
to real population data. If we ignore the reality for a 
distribution of population, the MSQ network includes a 
Sicrpinski gasket obtained by a special selection when 
each triangle, except the central one, is hierarchically di- 
vided. 






(a) triangle 



(b) square 
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(c) triomino chair tiling 
FIG. 1: Basic process of the division. 

The state-of-the-art geometric growing network mod- 
els [19-25] are summarized in Table I. The basic process 
for network generation is based on the division of a trian- 
gle or on the extension of an edge with a bypass route as 
shown in Fig. 3. For these models, we can also consider 
a mixing of dense and sparse areas of nodes by selecting 
a triangle or an edge according to a population in the 
territory. Although the geographical position of a new 
node has not been so far exactly defined in the geomet- 
ric processes, it is obviously different from that in the 
MSQ network. Moreover, in the Pseudofractal SF net- 
work [23], the link length tends to be longer than that 
in the MSQ network, because a node is set freely at an 
exterior point. Whereas the newly added three nodes ap- 
proach each other in the Sierpinski network [19] as the 
iteration progresses, they are degenerated (shrunken) to 
one node as in the random Apollonian network [20, 21]. 
Thus, we focus on the Apollonian network constructed by 
a biased selection of a triangle according to a population. 
In Sec. Ill, we compare the topological and the routing 
properties in the BA-like and the Apollonian networks 
with those in the MSQ network. 

The above geometric models generate SF networks 
whose degree distribution follows a power-law. It is 
known that a SF network is extremely vulnerable against 
intentional attacks on hubs [3, 4], in particular the tol- 
erance of connectivity is more weakened by geographical 
constraints on linking [30]. However, the MSQ network 
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FIG. 2: (Color online) Example of the proposed geographical 
network. The mesh data of population statistics (consisting 
of 8 2 x 10 2 x 4 = 25600 blocks for 80km 2 ) provided by the 
Japan Statistical Association are mapped onto each triangle 
by counting the number of people in the space. From light 
(orange) to dark (red), the gradation on a mesh block is pro- 
portionally assigned to the population. In the middle of the 
main island of Japan, almost all of the blocks consist of flat 
lands and mountains, and do not come in contact with the 
sea except in the left bottom curved white area. 



without high degree nodes has a quite different prop- 
erty. The MSQ network consists of trimodal low de- 
grees: k\ — 2, &2 = 4, and k% = 6 for an initial triangle 
(&2 = 3, and k^ = 4 for an initial square) configuration. 
Because of the trimodal low degrees without hubs, the 
robustness against both random failures and intentional 
attacks is maintained [28] at a similar level as the opti- 
mal bimodal networks [31] with a larger maximum degree 
k>2 = 0(\f~N) in a class of multimodal networks, which 
include a SF network at the maximum modality — > oo as 
the worst case for the robustness. 

In the MSQ network, the construction method defined 
by recursive subdivision of equilateral triangles or squares 
can be extensively applied to self-similar tilings based on 
polyomino [32], polyiamond [33], and polyform [34], as 
shown in Fig. 1(c). In the general construction on a 
polygon, some links are removed from the primitive tiling 
which consists of equilateral triangles or squares in order 
to be a specially shaped face such as a "sphinx" or "chair" 
at each time step of the subdivision. It remains to be 
determined whether or not the robustness of connectivity 
is weaker than that in our geographical networks based 
on equilateral triangles or squares. However, the path 
length becomes larger at least, as mentioned in Sec. Ill 
A. 



TABLE I: Geometically constructed network models. 





(a) Random Apollonian Net 



(b) Random Pseurdofractal SF Net 




a H 

(c) From a random Sierpinski gasket to the correspondig Sierpinski Net 

FIG. 3: Other geometric networks. From an initial triangle 
ABC, the following processes are repeated, (a) At each time 
step, a triangle is randomly chosen. Then, a new node is 
added (e.g., at the barycenter as a well-balanced position) 
and connected to the three nodes of the chosen triangle for 
the division, (b) Instead of a triangle, an edge is randomly 
chosen. Then, a new node is added at an exterior point, and 
connected to the two ends of the chosen edge, (c) After a 
random selection of a triangle, it is divided into small ones, 
and the center (shaded part) is removed to construct a random 
Sierpinski gasket. Then, the corresponding random Sierpinski 
network is obtained from the mapping of nodes and contacts 
to links between them. Note that the Sierpinski network has 
a similar structure to the Apollonian network rather than the 
Sierpinski gasket. 



III. ROUTING PROPERTIES 

A. Bounded shortest path 

The proposed MSQ network becomes the t-spanner 
[35], as a good graph property known in computer sci- 
ence: the length of the shortest distance path between 
any nodes u and v is bounded by t times the direct Eu- 
clidean distance d uv . A sketch of the proof is shown in 



the Appendix. Here, t is called stretch factor which is 
defined by a ratio of the path length (a sum of the link 
lengths on the path) to d uv . Figure 4 shows typical cases 
of the maximum stretch factor t = 2 in the MSQ network 
[28]. When the unit length is defined by an edge of the 
biggest equilateral triangle (or square), the path length 
denoted by a dashed line in Fig. 4 is 1 (or 2) while the 
direct Euclidean distance between the two nodes is 1/2 
(or 1). 



*) 




FIG. 4: Cases of the maximum stretch factor t — 2 for 
the shortest distance path denoted by a dashed line and 
the straight line directly passing the gap between two nodes 
marked by circles. 



More generally, if we construct a network by recur- 
sive subdivision to a self-similar tiling of a polyform [34] , 
e.g., polydraftcr: consisting of right triangles, polyabolo: 
consisting of isosceles right triangles, or domino: consist- 
ing of rectangles, then the stretch factor can be greater 
than 2 [see the U-shaped path in the right of Fig. 1(c)]. 
Thus, our network model based on equilateral trian- 
gles or squares is better for realizing a short routing 
path because of its isotropic property. In other ge- 
ometric graphs, the maximum stretch factor becomes 
larger: t — 27r/[3cos(7r/6)] sa 2.42 for Delaunay tri- 
angulations [36], and t = 2a > 4^3/3 « 2.3094 
for two-dimensional triangulations with an aspect ra- 
tio of hypotenuse/height less than a [37], whose lower 
bound is given for the fattest equilateral triangles. Al- 
though 0-graphs [38] with K non-overlapping cones have 
t = 1 / [cos(2ir / K) — sin.(2ir/K)] —J- 1 asymptotically as 
K — > oo, a large amount of O(KN) links is necessary, 
and some links may be crossed. In general graphs, even 
the existence of a bounded stretch factor is uncertain. 
On the other hand, in a SF network, the efficient routing 
[39] based on the passing through hubs has a stretch less 
than 2, which is defined by the ratio of the number of 
hops on the routing path to that on the shortest path. It 
can be implemented by a decentralized algorithm within 
small memory requirements. 

Since the geographical MSQ network is planar, which 
is also suitable for avoiding the interference among wire- 
less beams, we can apply an efficient routing algorithm 
[29] using only local information of the positions of nodes: 
neighbors of a current node, the source, and the terminal 
on a path. The online version has been developed in a 
distributive manner, in which necessary information for 
the routing is gathered through an exploration within a 
constant memory. As shown in Fig. 5, by using the face 
routing, the shortest distance path can be found on the 
edges of the faces that intersect the straight line between 
the source and the terminal nodes. Note that, in other 
decentralized routings without global information such 
as a routing table, some of them in early work lead to 
the failure of guaranteed delivery [40]; e.g., in the flood- 
ing algorithm, multiple redundant copies of a message 
are sent and cause network congestion, while greedy and 
compass routings may occasionally fall into infinite loops. 
On a planar graph, the face routing has advantages for 
the guarantee of a delivery and for the efficient search on 
a short path without the flooding. 
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FIG. 5: The shortest path obtained by the face routing on 
the edges of the shaded triangles that intersect the straight 
line between the source and the terminal nodes. Two paths 
of dashed segments in the opposite directions (for nl-nl' and 
v3-v3') are canceled, therefore the route does not pass the 
points vi and V3, and through only the other nodes (marked 
by open circles) on the triangles. Exceptionally, it takes the 
shortcut v5'-v6 to avoid the detour passing through V5 and 
V B ". 



B. Comparison of P(k) and P(Lij) 
networks 



with other 



We investigate the distributions of degree and of link 
length related to communication costs. In the following 
BA-likc networks, the positions of nodes arc fixed as same 
as in the original MSQ network, and only the connections 
are different. From the set of nodes in the MSQ network, 
a node is randomly selected as the new node at each 
time step in the growing process. On the other hand, the 
positions of nodes in the Apollonian network are different 
from the original ones because of the intrinsic geometric 
construction. 

Let us consider an extension of the BA model with 
both effects of distance [5-9] and population on linking. 
Until a size -/V is reached, at each time step, m links are 
created from a new node i to already existing nodes j. 
The attachment probability is given by 



IF, oc d tJ a pop^k], 



(1) 



where a, /3, 7 > are parameters to control the topol- 
ogy [5, 12], kj denotes the degree of node j, and d%j 
denotes the Euclidean distance between nodes i and j. 
The newly introduced term popj is not constant but may 



vary through time. If the nearest node from a point on 
the geographical space is changed by adding a new node, 
then the assigned population popj in the territory of the 
affected node j is updated. We set the average degree 
(k) = 2m = 4 that is the closest integer to (k) in the 
MSQ network. 

Each term in the right-hand side of Eq. (1) contributes 
to making a different topology, as the value of a, j3, or 7 is 
larger in the competitive attachments. As shown in Fig. 
6, a proximity graph is obtained in the BA-like:300,330 
networks by the effect of distance for a = 3, and hubs 
emerge near large cities in the BA-like:030,033,333 net- 
works by the effect of population for j3 = 3, while in the 
BA-likc:303,003 networks, a few hubs emerge at the posi- 
tions of randomly selected nodes in the early stage of net- 
work generation without any relation to the population. 
In the following discussion, we focus on BA-like:300,033 
networks, since they are typical with the minimum and 
the maximum degrees or link lengths, respectively, in the 
combination of a/?7 except 000,003 without both effects 
of distance and population. Figure 7(a) shows that the 
degree distributions P(k) follow an exponential decay in 
the BA-like:300 network, and a power-law like behavior 
in the BA-likc:033 and the Apollonian networks, which 
are denoted by a (red) solid line, (green) light dashed line, 
and (blue) dark dotted line, respectively. Figure 7(b) 
shows the distributions P(Lij) of link length counted as 
histograms in the interval 0.05. The Ly is normalized 
by the maximum length on the outer square or trian- 
gle. Some long-range links are remarkable in the BA- 
like networks, while they are rare in the quickly decaying 
distribution in the MSQ network. The average lengths 
shown in Table II are around 0.2 ~ 0.3 in many cases of 
the BA-like networks, and 0.02 ~ 0.14 in the Apollonian 
network. The original MSQ network has the smallest link 
length in less than one digit. As the size TV increases, the 
link lengths tend to be shorter due to a finer subdivision 
in all of the networks. In any case, the BA-like networks 
have longer links even with neighboring connections than 
the corresponding MSQ network, while the Apollonian 
network shows the intermediate result. Thus, the MSQ 
network is better than the Apollonian and the BA-like 
networks in term of the link length related to a commu- 
nication cost. 



C. Heavy- loaded nodes and links 

In a realistic situation, packets arc usually more often 
generated and received at a node, as the corresponding 
population is larger in the territory of the node. Con- 
sequently, the spatial distribution of the source or the 
terminal node is not uniformly random. Thus, we demon- 
strate how the traffic load is localized in the case when 
the number of generated and received packets at a node 
varies depending on a population assigned to the node. 

The traffic loads at a node i and through a link I are 
measured by the effective betweenness centralitics Bi and 
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FIG. 6: (Color online) Visualizations of the original MSQ, the 
Apollonian, and the typical BA-like networks for N = 100. As 
a common feature, many nodes concentrate on dark cloudy 
(red) areas with large populations. The positions of nodes 
are same in the MSQ and the BA-like networks, while they 
are slightly different in the Apollonian network, (b) Remark- 
able hubs exist at the center and four points on the outer 
lines, (c) Almost nearest connections are constructed by the 
attachment of distance, (d) Hubs emerge near large cities in 
the dark cloudy (red) areas: Kobe, Osaka, and Kyoto by the 
attachments of population and degree. These results with the 
nearest connections and hubs are similarly obtained in other 
combinations of afi-y for BA-like networks. Here, or 3 in the 
triplet denotes the parameter values for a, (3, and 7 in Eq. 

(I)- 



B[ , which are defined as follows [41]: 
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where b J k is the number of shortest distance paths between 
the source k and the terminal j, tf.ft) is the number of 
the paths passing through node i, and bj,(l) is the num- 
ber of the paths passing through link I. The first terms 
on the right-hand side of both Eqs. (2) and (3) are nor- 
malization factors. Although the measured node i and 
link I are usually excluded from the sum in the definition 
of betweenness centralities [42] , we include them tanking 
into account the processes for the generation and the re- 
moval of a packet in these measures Bi and B\ in order 
to investigate all of the traffic loads. 

Figures 8 and 9 show the link load measured by Bi for 
the following two selection patterns of packet generations 
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FIG. 7: (Color online) Distributions of (a) degree and (b) 
link length in the averaged 50 realizations of the Apollonian 
and the typical BA-like networks for N = 1000. (a) With 
power-law like behavior, there exists a few huge hubs around 
k = 10 2 ~ 10 3 at the vertical (green) light dashed bars for 
the BA-like:033 network, and ordinary hubs on the horizontal 
(blue) dark dotted line for the Apollonian network. In spite of 
the stochastic process, the locations of hubs are fixed e.g., at 
the center and four points on the outer lines in Fig. 6(b). Thus, 
the frequencies are high at some special degrees, (b) The 
histograms are shown with shifts for better discrimination. 
The gray dashed lines with filled triangles show P(Lij) in the 
original MSQ network. 



or removals. In the first pattern denoted by Pop, as the 
source or the terminal, a node is chosen with a probability 
proportional to the assigned population to the node (left 
figures), while in the second pattern denoted by Rand, it 
is chosen uniformly at random (right figures). Therefore, 
the sum in Eqs. (2) or (3) has a biased frequency for each 
pair of nodes by their populations in Pop, while the sum 
corresponds to the simple combination of nodes in Rand. 
In Pop, heavy-loaded links denoted by thick (blue) lines 
between the cloudy (red) areas with large populations are 
observed: the routes between well-known cities Osaka 
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TABLE II: Average link length (Ly) in the original MSQ, the 
Apollonian, and the BA-like networks. The (Lij) becomes 
longer in the order: the MSQ, the BA-like:300 or the Apollo- 
nian, and the BA-like:033 networks. This order is consistent 
with the positions of peaks and the widths of P(Lij) in Fig. 
7(b). Note that BA-like:000,003 networks are not geographi- 
cal, but exceptional, because they have no effects of distance 
and population on linking. Here, the triplet of or 3 in the 
first column denotes the values of a, ft, and 7. 



and Kyoto (in the enclosed bold lines) are remarkable 
with thicker (blue) links than in Rand. Note that the 
number N — 100 of nodes is due to a comparatively clear 
appearance of the difference in two patterns. At a larger 
N, the packet transfer is already biased by the spatial 
concentration of nodes, and the assigned population in 
each territory of node is well-balanced by the division 
process. Then, the difference of heavy-loaded links in 
Rand and Pop tends to disappear. However, it is worth 
considering such biased selections of packets in a realistic 
problem setting, since the spatially localized positions of 
the heavy-loaded links are not trivially predictable from 
the results for the usually assumed Rand. 

We further investigate the maximum betweenness in 
the scaling relation N s which affects traffic congestion 
[43]. Here, the betweenness is defined by the number of 
passings through a link or a node on the shortest distance 
paths. In general, a smaller S yields a better performance 
for avoiding traffic congestion. The exponent value S de- 
pends on both a network topology and a routing scheme. 
In this paper, assuming the shortest distance paths ob- 
tained by the face routing in Sec.IIIA, we focus on an 
effect of the topologies on the scaling relations in the 
MSQ, the Apollonian, and BA-like networks. 

Figure 10 shows the scaling relations of the maximum 
node betweenness. They are separated into two groups 
of S = 1.67 - 1.77 and 6 = 1.87 - 1.97. The MSQ 
network with only trimodal low degrees is located on the 
baseline. The order of thick lines from the bottom to the 
top fairly corresponds to the increasing order of largest 
degrees in these networks, as shown in Fig. 7(a), be- 
cause more packets tend to concentrate on large degree 
nodes as the degrees are larger. Figure 11 shows the 
scaling relations of the maximum link betweenness. The 
thick lines from the top to the bottom in the inverse or- 
der to that in Fig. 10 show that packets are distributed 
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FIG. 8: (Color online) Visualizations of the link betweenness 
centrality Bi in a MSQ network for N - 100. Pop (left) and 
Rand (right) denote the selection patterns of a source or a 
terminal node, which are chosen proportionally to the popu- 
lation in the territory of a node, and uniformly at random, 
respectively. From the thin (green) to the thick (blue) verti- 
cal or horizontal line, the gradation is proportionally assigned 
to the value of Bi. The enclosed bold line in Pop emphasizes 
the parts of large Bi with heavy load, which are remarkably 
shown as thick (blue) links on the paths connected to dark 
cloudy (red) areas with large populations in a diagonal direc- 
tion. The enclosed bold lines in Rand emphasize the thick 
(blue) lines on light gray (orange and white) background for 
underpopulated areas, as the other differences between the 
left and the right figures. Thus, in Rand, some thick (blue) 
lines for large Bi have no relation to the population. 
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FIG. 9: (Color online) Visualized examples of the link be- 
tweenness centrality Bi in the BA-like:300 (top) and the Apol- 
lonian (bottom) networks for N — 100. The enclosed bold 
line in Pop emphasizes that the thick (blue) links for large Bi 
are on the paths connected to dark cloudy (red) areas with 
large populations. Pop (left) and Rand (right) denote the 
two selection patterns of a source or a terminal node which is 
chosen proportionally to the population assigned to a node, 
and uniformly at random, respectively. 



on many links connected to large degree nodes as the 
degrees are larger. However, all lines lie almost around 
the intermediate slopes of S = 1.6, especially for a large 
size N, the maximum link load is at a similar level in 
the lowest case of the Apollonian network and the high- 
test case of the BA-like:300 or the MSQ network. These 
results in Figs. 10 and 11 are also consistent with the 
increasing orders of the maximum Bi and Bi in Table 
III. Note that there exists no remarkable difference be- 
tween the two patterns of Rand (marked by triangles) 
and Pop (marked by inverted triangles) for both scaling 
relations of node and link loads except in the Apollonian 
network (dashed-dotted line). In summary, the MSQ net- 
work yields a better performance with a lower maximum 
load than the other geographical networks, therefore it is 
more suitable for avoiding traffic congestion. 
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0.854 
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0.362 


0.392 


0.134 


0.180 


033 


0.657 


0.710 


0.056 


0.099 


300 


0.282 


0.235 


0.121 


0.121 


303 


0.444 


0.532 


0.073 


0.107 


330 


0.397 


0.337 


0.121 


0.143 


333 


0.634 


0.620 


0.106 


0.072 


Apollonian 


0.295 


0.278 


0.056 


0.059 


MSQ 


0.227 


0.259 


0.137 


0.185 



TABLE III: Maximum betweenness centralities in the BA- 
like, the Apollonian, and the MSQ networks for N = 100. 
The triplet of or 3 in the first column denotes the values of 
a/?7 in the BA-like networks. We note that the MSQ, the BA- 
like:300, the Apollonian, and the BA-like:033 networks are in 
increasing order of max Bi and decreasing order of max Bi . 



IV. CONCLUSION 

Both mobile communication and wide-area wireless 
technologies (or high-speed mass mobilities) are impor- 
tant more and more to sustain our socio-economic ac- 
tivities, however it is an issue to construct an efficient 
and robust network on the dynamic environment which 
depends on realistic communication (or transportation) 
requests. In order to design the future ad hoc networks, 
we have considered geographical network constructions, 
in which the spatial distribution of nodes is naturally 



determined according to a population in a self-organized 
manner. In particular, the proposed MSQ network model 
is constructed by a self-similar tiling for load balancing of 
communication requests in the territories of nodes. On a 
combination of complex network science and computer 
science approaches, this model has several advantages 
[28]: the robustness of connectivity, the bounded short 
paths, and the efficient decentralized routing on a planar 
network. Moreover, we have numerically shown that the 
MSQ network is better in term of shorter link lengths 
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FIG. 10: Scaling relations TV* of the maximum node between- 
ness in the MSQ (solid line), BA-like:300 (dotted line), the 
Apollonian (dashed-dotted line), and the BA-like:033 (dashed 
line) networks. The triangle and inverted triangle marks cor- 
respond to Pop and Rand, respectively. The thin lines with- 
out any marks guide the slopes 5 — 1.67, 1.77, 1.87, 1.97 from 
the bottom to the top. 
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FIG. 11: Scaling relations N of the maximum link between- 
ness in the MSQ, the BA-like, and the Apollonian networks. 
The thick lines and marks are the same as in Fig. 10. The 
thin lines without any marks guide the slopes S — 1.5, 1.6, 1.7 
from the bottom to the top. 



than the geographical Apollonian and the BA-like net- 
works which have various topologies ranging from river 
to SF networks. As regards the traffic properties in the 
MSQ network, the node load (defined by the maximum 
node betweenness) is lower with a smaller 5 in the scaling 
relation N s , although the link load is at a similar level as 
that of the other geographical networks. Therefore, the 
MSQ network is more tolerant to traffic congestion than 
the state-of-the-art geographical network models. 

In a realistic situation, packets arc usually more often 
generated and received at a node whose population is 



large in the territory. Thus, we take into account spa- 
tially inhomogeneous generations and removals of pack- 
ets according to a population. Concerning the effect of 
population on the traffic load, the heavy-loaded routes 
with much throughput of packets are observed between 
large population areas especially for a small size N. The 
heavy-loaded parts are not trivial but localized depend- 
ing on the geographical connection structure, whichever 
a spatial distribution of nodes or biased selections of the 
source and the terminal is more dominant. In further 
studies, such spatially biased selections of packets should 
be investigated more to predict overloaded parts on a re- 
alistic traffic, since the selections may depend on other 
economic and social activities in tradings or community 
relationships beyond the activities related only to pop- 
ulation, and probably affect the optimal topologies for 
the traffic [11, 13, 15]. The optimal routing [43] instead 
of our shortest distance routing is also important for re- 
ducing the maximum betweenness (Bi or B{) and the 
exponent 5 in the scaling relation. In the optimizations 
that include other criteria [9, 14], we will investigate the 
performance of the MSQ networks. Related discussions 
to urban street networks [44-46] are attractive for inves- 
tigating a common property which has arisen from the 
geographically pseudofractal structures. 
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Appendix 

We provide a sketch of the proof for the i-spanner prop- 
erty in the MSQ networks. The following is similarly dis- 
cussed for other polygons by exchanging the term "trian- 
gle" or "triangulation" with other words (e.g., "square" 
or "polygonal subdivision"), except that the maximum 
stretch factor t may be greater than two for a general 
polygon. Remember that the self-similar tiling is ob- 
tained from the division by a contraction map as shown 
in Fig. 1. 

Let T be a given triangulation and L st be a line seg- 
ment interconnecting two vertices so and te in T. It inter- 
sects many triangle faces in T. Let (so = vo, v\, . . . , v n = 
te) be those ordered vertices. Now, we define a path us- 
ing triangular edges. For any two consecutive vertices Vi 
and Vj-|_i there is a unique triangle which contain both of 
them. In other words, Vi and Vi+\ are the entrance and 
the exit of the line segment to the triangle. For the pair 
(vi,Vi + i) there are two paths, clockwise and anticlock- 
wise, on the triangle. We take the shorter one. In this 
way, we can define a path on the given triangulation. If 



the path contains duplicated segments with the opposite 
directions as shown in Fig. 5, we remove them. Then, 
we have a path such that each edge of the path is some 
triangular side. When the two consecutive edges on the 
path belong to the links of a same triangle, we take a 
shortcut by directly passing another edge (see two edges 
v5'-v5" and v5"-v6 are replaced with vb'-v 6 on the path 



in Fig. 5). Since T is planar and there is no node inside 
each triangle, the length of this path is the closest to 
that of the line segment L st , therefore the shortest in all 
other routes. As is easily seen, in each triangle, the path 
length between Vi and Vi+\ is at most twice longer than 
the Euclidean distance between them. The worst case of 
the maximum stretch factor t = 2 is illustrated in Fig. 4. 
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